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Abstract 

Photonic SUSY two-loop corrections to the muon magnetic moment are 
contributions from diagrams where an additional photon loop is attached 
to a SUSY one-loop diagram. These photonic corrections are evaluated ex- 
actly, extending a leading-log calculation by Degrassi and Giudice. Com- 
pact analytical expressions are provided and the numerical behaviour is 
discussed. The photonic corrections reduce the SUSY one- loop result by 
7 . . . 9%. The new terms are typically around ten times smaller than the 
leading logarithms, but they can be larger and have either sign in cases 
with large SUSY mass splittings. We also provide details on renormaliza- 
tion and regularization and on how to incorporate the photonic corrections 
into a full SUSY two-loop calculation. 
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1 Introduction 



The muon magnetic dipole moment belongs to the most precisely known observ- 
ables in particle physics. The anomalous magnetic moment = {g — 2Jl/2 has 
been determined by an impressive series of measurements at BNL to [ij^l 

a^fP = (11659 208.9±6.3) X 10"^°. (1) 

At this level of precision, is sensitive to all interactions of the Standard Model 
(SM) as well as to hypothetical new particles at the electroweak scale. 

In recent years, the precision of the Standard Model (SM) theory evaluation 
has reached a similar level. The hadronic vacuum polarization contributions can 
be related to the cross section for e^e^ — )■ hadrons in a theoretically clean way, 
and crucial recent progress has been achieved on the experimental determination 
of this cross section by the SND, CMD-2, KLOE, and Babar experiments [3H6]. 
The quality of the current e^e~ — )■ hadrons data is such that all recent theory 
evaluations of the corresponding contributions to agree fairly well, see e.g. 
[7Hl2]l!l 

The hadronic light-by-light contributions have been scrutinized by many 
groups in the past 15 years, and while an ever better understanding has been 
achieved, the central value of the result has remained relatively stable. Re- 
cently, three groups have joined forces and published a common value, aj^^ = 
10.5(2.6) X 10^^° [17J, where the errors have been enlarged in order to cover the 
results obtained from the different approaches. Compatible results have been 
obtained even more recently in Refs. [T0|[T8]. 

As a result of the recent progress, the SM theory prediction for now has 
an even smaller error than a^^P, calling for a new experiment. For reference we 
use the value from [12], 

aj^ = (11 659 183.4 ±4.9) X 10"^°. (2) 

The deviation from the experimental value is 

Aa^(exp - SM) = (25.5 ± 8.0) x 10"^° . (3) 

All the other mentioned recent e^e^-based evaluations lead to similar deviations 
in the range (25.5 . . . 31.6) x 10~^° with combined errors in the range (7.9 . . . 9.0) x 
10-1°. 

^The change in the number compared to Ref. [I] is due to a new PDG value for the magnetic 
moment ratio of the muon to proton [5] . 

^In principle, part of the e^e~ — > hadrons cross section could be obtained in an alternative 
way from hadronic r decays |13j . However, the required isospin breaking effects are difficult 
to quantify at the currently needed level of precision [10 11141116] . Hence, in order to obtain 
the most reliable results, most analyses do not use r decays. Nevertheless, the most recent 
analysis employing an improved understanding of isospin breaking effects shows a marginal 
consistency between the e"'"e~-based and the r-decay based evaluations of a^. 



2 



This 3-4(7 deviation constitutes a tantalizing hint for physics beyond the SM. 
Although the deviation is almost twice as large as the SM weak contributions aj^^'^^ 
and contributions from hypothetical heavy particles with mass M are typically 
suppressed oc {Mw / M)^ aj^'^'^^ , there is a variety of models that could explain 
it [in]. Supersymmetry (SUSY) is a particularly promising example (for a review 
see 120]), owing to an enhancement by tan/3, the ratio of the two Higgs vacuum 
expectation values. 

Even more importantly, the precision of the deviation ([3]) implies significant 
constraints on the parameters of many new physics models. In supersymmetry, 

permits to derive mass limits can which cannot be obtained from other ob- 
servables [21]; is also one of the most important quantities in recent global 
analyses of the parameter space of various supersymmetric models [221425] . Even 
in the LHC era, will remain a highly useful complementary observable. It 
will provide a benchmark for models, help selecting between different SUSY sce- 
narios, eliminate ambiguities, and it will improve parameter determinations. For 
instance, combining LHC data with will significantly reduce the tan /3 uncer- 
tainty [261128]. 

Given the current situation and the usefulness of a^, every effort should be 
made in order to improve the experimental and theoretical precision of a^. In the 
near future, the precision of Aa^(exp — SM) will indeed further increase. New 
analyses of e~^e^ — t- hadrons data are in the pipeline, and a new CMD-3 experi- 
ment in Novosibirsk is planned. All of these will directly feed into a more precise 
evaluation of the hadronic vacuum polarization contributions. Furthermore, a 
new, improved measurement of itself could be carried out at Fermilab [29] 
or possibly at J-PARC [30]. Building on the BNL experiment, whose final un- 
certainty is still statistics dominated, the Fermilab measurement could reach a 
reduction of the uncertainty down to 1.6 x 10^^'^. 

With these expected improvements, the theory error from contributions from 
physics beyond the SM, in particular from supersymmetry, becomes more promi- 
nent. The current theory error of the SUSY contributions has been estimated 
to 3 X 10^^° [20], which is larger than the precision goal of the Fermilab 
measurement. 

The one-loop contributions from supersymmetric particles have been com- 
puted a long time ago [3T1[32], and the SUSY two-loop corrections to SM one- 
loop diagrams (e.g. SM-diagrams with insertions of closed chargino/neutralino 
or stop/sbottom loops) are completely known [S5HS5] . However, the two- loop 
corrections to SUSY one-loop diagrams are not known, only two dominant parts 
have been identified: large QED-logarithms [36] and (tan /3) ^-enhanced correc- 
tions |37j. A subclass of the remaining diagrams has been considered in Ref. [38j . 
The unknown two-loop contributions lead to the theory error mentioned above. 

In the present paper we consider the photonic two-loop contributions, defined 
as contributions from diagrams where a photon loop is attached to a SUSY one- 
loop diagram. In Ref. [36] , Degrassi and Giudice showed that these contributions 
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lead to large QED- logarithms of the form log(M/m^), where M is the mass scale 
of the new particles. They evaluated the logarithms using elegant effective field 
theory and renormalization group techniques. These logarithms reduce the one- 
loop contributions by 7 . . . 9% for M between 100 and 1000 GeV. It is clear that 
there is an ambiguity in the choice of M, in particular in the case of a rather split 
spectrum. In the present paper we evaluate the photonic two-loop corrections 
exactly. In this way, we resolve the ambiguity, and we obtain all logarithms 
of ratios of different heavy masses and the associated non-logarithmic terms. 
Our results are derived in a generic model which covers the case of the minimal 
supersymmetric SM (MSSM) but also of a wider class of models. We take into 
account issues such as the choice of dimensional regularization versus dimensional 
reduction, and we set up the calculation in such a way that our results will be 
useful building blocks for a full two- loop computation of a^^^^ ■ 

The outline of the paper is as follows. In section |2] we describe the setup 
of our calculation and provide the necessary one-loop contributions up to the 
order e on the regularized level. In section [3] we list the contributing two-loop 
diagrams, classify them according to their ultraviolet (UV) and infrared (IR) 
divergences, and explain our method to compute them. Section H] is devoted 
to renormalization and the analysis of the counterterm diagrams. We provide 
details on the cancellation of UV and IR divergences and the regularization- 
scheme dependence. In section [5] we provide our results as a compact analytical 
formula, and discuss several numerical examples. In section |6] we discuss how 
our results can be used as building blocks in a full MSSM two-loop calculation. 
Section [7] contains the conclusions. 



2 Generic model and one-loop contributions 

The genuine SUSY one-loop contributions to in SUSY extensions of the SM are 
given by the two kinds of diagrams in Fig. [TJ The loops involve either a chargino 

and a sneutrino i/^ or a neutralino and a smuon jl. At the two-loop level it 
is useful to classify the SUSY contributions into those from diagrams where the 
/i-lepton number is only carried by muon or muon-neutrino, and diagrams where 
the /i-lepton number is carried also by jl and/or z/^. The second class can be 
interpreted as the two-loop corrections to the SUSY one-loop diagrams in Fig. [1] 

In the present paper we aim for the computation of the photonic contributions 
of this second class, i.e. of two-loop diagrams obtained from attaching a photon 
loop in all possible ways to the SUSY one-loop diagrams. The relevant Lagrangian 
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Figure 1: The two one-loop diagrams involving SUSY particles. In this and all 
following figures, thin lines denote the photon and muon, thick lines denote the 
chargino/neutralino, thick dashed lines the sneutrino/smuon. 

is given by 

-^QED + ^mat + -^int, (4) 

+ lx^{^^ - m^o)x' + \d^i>\' - ml\9\\ (5) 
Ant = z>^^(c2Pl + CRPR)fi + fi^^inlPL - nRPji)fi + h.c. (6) 

Here £qed is the QED Lagrangian for photon and muon /i; x~ and x° are the 
chargino and neutralino spinor fields, D^^ is the covariant derivative in QED and 
Pl,r = ^(1 -F 7^)- This Lagrangian contains only photon and muon plus one 
generic chargino, neutralino, sneutrino, and smuon. For the following, it will be 
useful to adopt the viewpoint that CQED+c,n describes a simple generic model 
which extends QED by two charged and two uncharged particles. The couplings 
of the SUSY particles to the muon are given by free parameters c^'^ and n^'^, 
all other terms in the Lagrangian are determined by QED gauge invariance. The 
conventions of the couplings c, n are chosen in the same way as in Ref. [39]. In 
terms of this generic model, the one- loop contributions in Fig. [1] and the photonic 
two-loop corrections (plus renormalization) are defined as the contributions of 
C(c^,n^) and 0{ac^ ,an'^), respectively, with the fine structure constant a = 
eV47r. 

By specializing the couplings, the generic model can describe photonic cor- 
rections in underlying models such as the MSSM, but also in non-minimal SUSY 
models, and even in non-supersymmetric models, if they happen to contribute to 

by diagrams of the form in Fig. [H For instance, in the MSSM, the couplings 
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2^'^ and n^'^ have the following values: 

Cfc = -92Vki, (7a) 

cf = y^.Uk2, (7b) 

nt = -^(aiNa + g2N^2Wmi - y.N^.Uf'^,*, (7c) 

V2giNiiUf'^2 + 2/MiV,3[/'^„i, (7d) 



where the chargino, neutralino and smuon indices can take the values k G {1,2}, 
i E {1, . . . , 4}, m G {1, 2}, and where the notation of Ref. ^20j has been used for 
gauge and Yukawa couplings and mixing matrices. 

As a first step we present the result of the one-loop diagrams of Fig. [1] These 
diagrams will be needed later as counterterm diagrams, multiplied with 1/e di- 
vergences from renormalization constants, where D = 4 — 2e in dimensional 
regularization. Hence, we do not only need the finite part of the result but also 
the C(e)-part. 

This (9(e)-part can in principle depend on the regularization scheme. We 
consider two schemes: dimensional regularization (DREG) and dimensional re- 
duction (DRED). As discussed e.g. in |10], it is crucial to define the treatment 
of external vector bosons in dimensional schemes. Our definitions are as follows: 
In DREG, the vertex function T^p^p is computed with a D-dimensional photon, 
then the D-dimensional projector to extract defined in is applied. In 

DRED, the vertex function F^^^p is computed with a (quasi-) 4- dimensional pho- 
ton, then the photon is projected onto the D-dimensional subspace, Qp^T ^^^a" in 
the notation of |40j, and finally the same Z)-dimensional a^-projector is applied|§ 

With these definitions, it is clear that the regularized results for the diagrams 
in Fig. [T]are equal in DREG and DRED, even at 0{e). The results read 



^0 ^ 1 ml 
167r2 ml 



|{l^^^f(x) + ^S^^f(a:)}, (8a) 
{-l^^^f(x) + ^S^^,^(x)}, (8b) 



with the kinematic variables x = m'^±/ml , or x = m?o/m\, respectively, and 
the coupling combinations 

= + |n«p, = (9b) 



•^Note that in DRED, for a quasi-4-dimensional photon, the covariant decomposition of 
^fj.p.A'^ has more terms than in DREG, and that different prescriptions to project out can 
lead to different results at 0{e). 
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The e-dependent loop functions are decomposed as 

J-f (x) = - eLKj] + ei^?(x), (10) 

^/"(^) = F.r{x) [1 - 6 L(mJ)] + 6 i^f (x), (11) 

where we have used the abbreviation 
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171 

Um')=\og- (12) 

/^DREG 

with the dimensional-regularization scale /iDREG? and the well-known functions 

2 

Ff (x) = r- \2 + 3x- + + &X logxl , (13) 

(1 — xp 

F^{x) = ^^^^7^ [ - 3 + 4x - - 2 log x] , (14) 

F^(x) = - — [1 - 6x + 3x^ + 2x^ - 6x^ logxl , (15) 
(1 — x)* 

^2^^^^) = (i_^^)3 [1 - + 2x logx] , (16) 

normalized such that F-{1) = 1. The functions for the e-dependent parts are 
defined as 

F^Ax) = (x) ( -^^ + 6^^ + 15a: + 2-6xlogx\ ^ x^-8x-4^ ^^^^ 

= (X) (z2£l±^^±l^li2I£) + (18) 

= F,"(.) ?£!±±L+_Lzl£i^^ _ !£±1, (20) 

and are normalized to F-^l) = 0. The (9(e°)-part of this result has been presented 
in a form similar to the one in [39] . 

Our definition of the coupling combinations A, B reflects an important phys- 
ical property. The muon anomalous magnetic moment, like the muon mass, 
corresponds to a chirality-flipping operator, and all contributions must be pro- 
portional to a chiral-symmetry violating parameter. In the SM, the MSSM and 
many other models, chiral symmetry in the muon sector is broken only by the 
muon mass and its Yukawa coupling. In the remainder of the article we assume 
nothing about the model underlying the Lagrangian (jl]) except that it has this 
same property. 
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With this assumption, all contributions to are of the order m'^^/M^jj^y in 
the muon mass, like in the MSSM p ^ [20 | l32]. Furthermore, the combinations 
A'"'^ and B'"'^ must have the behaviour 

j^C,N^ l^CN ^ ^Q^g^^ ^ (^(^2) ^21) 

in the muon mass. In other words, either or must be proportional to 
the muon mass, and similar for n^, n^. It is noteworthy that an equivalent 
assumption has been made implicitly in Ref. in the computation of the 2-loop 
QED-logarithms in new physics models using effective field theory techniques. In 
the matching between the full and the effective theory Ref. [36] assumed that one 
can pull a factor out of the new physics contributions to a^. If that were not 
the case, additional large QED-logarithms would appear in the matching. 

For the case of the MSSM, Eq. shows that the couplings indeed have the 
assumed behaviour, for all values of the chargino, neutralino and smuon indice^. 
For phenomenology it is important that i?*"'^ are enhanced by a factor tan /3. In 
a large part of the MSSM parameter space, the i3'-^-term in constitutes the 
dominant SUSY contribution to a... 



3 Two-loop contributions 

We are interested in the photonic two-loop corrections to Fig. [H i.e. the two-loop 
contributions of 0{ac^, an"^) in our generic model. We denote the corresponding 
genuine two-loop diagrams by af^'^'^^^l They are shown in Figs. Ej [3l divided 
into diagrams with chargino/sneutrino or neutralino/smuon exchange. All dia- 
grams are overall UV-convergent, but a useful classification according to divergent 
subdiagrams is possible. 

• The diagrams in Figs. [2^ and[3K contain a UV-divergent SUSY subdiagram, 
inserted into a QED diagram. 

• The diagrams in Figs. [2b and [3b contain a UV-divergent QED subdiagram, 
inserted into a SUSY diagram. There are more neutralino/smuon diagrams 
because of the presence of vertices involving two photons. 

• The diagrams in Figs. [2)3 and[3|3 contain no UV-divergent subdiagrams and 
are finite. 

There are also IR divergences in the third diagrams in Figs. ^ and [3|i, which 
are the only diagrams where both ends of the photon propagator are attached to 
the external muons and with a subdiagram that contributes to a^. 

"^In order to see this, note that the niuon Yukawa couphng and, for each to, one of the 
mixing matrix elements C/'^mi or U^m2 are proportional to to^. 
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Figure 2: Photonic two-loop corrections to chargino one-loop diagrams. In this 
and all following figures, it is understood that many graphs (here all graphs except 
the 3rd and 6th) appear twice, with the SUSY loop either on the upper or on the 
lower muon line. These graphs are displayed only once. 



(a) 




Figure 3: Photonic two-loop corrections to neutralino one-loop diagrams. 
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Figure 4: The QED counterterm diagrams and relevant renormalization con- 
stants. In the first and third diagram also 5e^^ ' ' would appear, but 

these vanish. 

We have evaluated the diagrams using DREG gulator for UV and IR 

divergences. In the next section we will discuss how the UV and IR divergences 
of the different classes of diagrams are cancelled by renormalization, and we will 
discuss the transition to DRED. Technically, we have employed computer algebra 
programs based on the packages FeynArts |12] (using both a custom model file 
and the standard MSSM model file [43j) and TwoCalc [44j (with a custom routine 
for the integral reduction). 

We have applied a large mass expansion [15] to all diagrams, where the muon 
mass is treated as small and all other masses as large. This results in a separation 
of heavy and light scales, and after reduction to master integrals only two kinds 
of loop integrals can appear: 

• products of a light one-loop integral (depending only on m^) and a heavy 
one-loop integral (independent of nif^): these are the integrals in which the 
large log m^-factors considered in [36] but also other terms are generated. 

• heavy two- loop integrals (independent of m^, at most multiplied with poly- 
nomials in m^): no large logarithms are generated here; these contributions 
are new. The heavy two-loop integrals have no external momentum and 
can depend only on two different mass scales. 

The structure of these integrals makes clear that log m^-factors appear in all 
diagrams in which the internal photon couples to at least one muon. Furthermore, 
the diagrams of Figs. [2^ and [3^ contain terms of order (m^/MsusY)°, while all 
other diagrams are suppressed by the small ratio {m^/Msusy)'^- The final result 
of all diagrams can be expressed as a rational function of the particle masses, 
and logarithms and dilogarithms. Below we will present the result obtained after 
renormalization. 

4 Renormalization and counterterms 

Apart from the actual two-loop diagrams, two kinds of counterterm diagrams 
can contribute at 0{an'^, ac^), and we denote the corresponding contributions as 
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Figure 5: The non-vanishing SUSY counterterm diagrams and relevant renormal- 
ization constants. 

^ct,QED _|_ ^ct,susY_ rjj^g QED counterterm diagrams in Fig. HI a^*'^^°, arise from 
renormahzation of the QED quantities in the model Lagrangian and involve the 
renormalization constants 



±,0 



(22) 



for muon mass and field renormalization, charge and photon field renormal- 
ization. The superscripts indicate that these QED renormalization constants 
need to be evaluated at (9(c^,?t,^), from SUSY one-loop diagrams involving 
a chargino-sneutrino or neutralino-smuon loop. This implies in particular 

(x±-") r (■.±fi^ 



5Z 



^g(x ' ) = 0. The values of the renormalization constants, and the 



physical meaning of the renormalized parameters, are fixed by the choice of a 
renormalization scheme. All QED quantities, in particular drrvj^ and SZ]^ 
must be defined in the on-shell scheme in order to guarantee the correct relation 
between the 3-point function and a^. 

The SUSY counterterm diagrams in Fig. O a^*'^^^"^, arise from renormaliza- 
tion of the SUSY quantities and involve the renormalization constants 



(23) 



for SUSY mass and coupling renormalization and muon field renormalization. 
Here all renormalization constants need to be evaluated at 0{a), i.e. from di- 
agrams with a photon loop. Mass and field renormalization constants for the 
neutral particles as well as ^e^'''-* and 5Z'^2^ would be zero and do not have 

to be included. Field renormalization constants 5Zf7,-, could be introduced in 
order to cancel UV divergences of individual two-loop diagrams. Since they drop 
out in the end, we ignore them here. 

Again, dZ^^"* has to be defined in the on-shell scheme. With this choice the 
photonic (i.e. 0{q)) contribution to the muon field renormalization constant is 
UV and IR divergent, and it can be decomposed as 



MM M ^yj- 



4 + 3L(mf 



(24) 



11 



The typical transcendental constants '~)e and log 4-7?, which cancel in the end, are 
ignored throughout, and L(m^) has been defined in Eq. (fT2l) . 

The choice of the renormalization scheme for the other parameters is more 
delicate. As discussed in Sec. |2] the photonic corrections appear as a subset 
of a larger class of two-loop contributions to in models such as the MSSM. 
In this context, the renormalization constants fl23|) are contributions to the full 
renormalization constants, 5m^"± = 5m^^ + ^ffi^'^^^^^'^'^^ etc. This split is not 
unique, and by choosing a scheme for the constants (!23l) we effectively define 
precisely what we mean by "photonic corrections" . 

A natural choice for the mass counterterms are the on-shell values of the 
photonic contributions to the self energies, so we choose 5rn:^± = T^^lij) = 'm^±) 

and 5m? = t!'^\p'^ = m?). The couplings c, n appear only via the combinations 
A'"''^ and m^B'"'^ , see Eq. ([9]). Their renormalization transformations c — )■ c + 
(5c*^'^^ and n — )■ n + 5n'^'~'^ can be equivalently written as 



In the class of models discussed in Sec. Ej the combinations A^''^ , B^'^ depend 
only on details of the underlying theory, while the muon mass is a low-energy 
quantity. Hence a natural choice is to renormalize A^'^ and B^'^ in the MS- 
scheme and the muon mass in the on-shell scheme. Implicitly, this choice defines 
a renormalization scheme for the couplings c, n. We call the scheme defined in 
this way the "on-shell muon mass scheme" 

As the major advantage, in this on-shell muon mass scheme the counterterm 
6{c^c^)'^'^^ / {c^c^) is set equal to drnf^'^ /m^ + dB'^^^'^B'" and similarly for n^n^, 
where the first term contains a large QED logarithm. This is important because 
in the considered class of models, c^c^ is proportional to the muon mass, and 
the full renormalization constants will satisfy 6{c^c^y^^^ /{c^c^) = 6771^^^/111^ + . . . 
and thus contain the same large QED logarithm. 

Rewriting the renormalization transformations in terms of differential oper- 
ators, we find a compact expression for the SUSY counterterm contributions in 

^In other words, in the on-shell muon mass scheme, the MS-conditions on A'-''^ and B^'^ 
imply that 5(|c^p + |c^p)^^^ is a pure MS-quantity while (5(c^c^)('''^ is not (and similarly for 
c — ^ n). Note that this is consistent at the considered order in m^, neglecting terms of ©(m^), 
but it makes difficult to individually solve for Sc^ and Sc^ as long as it is not specified which 
factor among c^, is proportional to m^i. 



(25) 
(26) 
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our on-shell muon mass scheme, 



^ct.SUSY 

V 



5Z^;\af+al') 



6m 



(7,fin) 



i=L,R 



r (7, fin) 



m„ 



(27) 



where it is understood that the partial derivatives d^i, d^i act on the couphngs 
within A^'^ and B^'^ . In spite of not using the MS-scheme for c, n we find it con- 
venient to express the counterterms in terms of the would-be MS-renormalization 
constants. For convenience, we provide the explicit values for the renormalization 
constants both in DREG and in DREE0 



5m 



(7) 

± 



ml 



c 



.L,R 



6m 



iiM) 
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a 
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a 
An 



+ 3L(m2±; 



- + SLfrnf 



3 
e 
3 
e 

3 
Ye 

3 



-A-e 



DRED 



-^DRED 



7DRED 



m,. 



An L 



m, 



4]- 



(28) 
(29) 
(30) 
(31) 
(32) 



The quantity 6'dred = in DREG and 6'dred = 1 in the case of DRED. 

We have checked that with these definitions all UV divergences cancel in 
the sum of the two-loop diagrams of Figs. [2^, [3^ and the QED counterterm 
diagrams with the appropriate counterterm insertions of (!?(n^,c^). Likewise, all 
UV divergences cancel in the sum of the two-loop diagrams of Figs. [Sb and the 
counterterms a^*'^^^^. Furthermore, we have checked that the UV divergences 
cancel for each diagram (plus corresponding counterterm diagram) separately, if 
X^- and /i-field renormalization is taken into account. 



See [46l|47] for a more general analysis of the transition from DREG to DRED. 
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The IR divergences cancel only in the following combinations: 




(34) 

The unsuppressed terms of order (m^/MsusY)° cancel between the two-loop 
diagrams of Figs. [21 [3] and the QED counterterm diagrams. 

Hence, after renormalization we obtain a UV and IR finite result of the order 
(m^/Msusy)^ for the photonic two-loop contributions to a^. 



5 Results 

The full result for the photonic two-loop corrections to a^, defined by the sum 

of the two-loop diagrams of Sec. El a^'^^'^^^'^\ and the counterterm contributions 
^ct,QED _|_ ^ct.susY^ ^^^^ ^ quite compact analytical form: 



,X± (7) 



1 a 



,x° (7) 



1 a 
IGtt^ 4:71 ml 



2m^± 



12 ^ ^ ' 3 



m, 



B^F^ix) 1 16 log ^ 



47 .c^r, N 122mv±^^^^, , 



-( ^A^ F^ {x) + 2m^±B^ F!^ {x) )L(m^ 



C TpCl 



- ^^^i^i^(^) + "^B^F.^x) ] 16 log 



9 



+ (-^^Vr(x))LK) 



(35) 



(36) 
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where the kinematic variables are again defined as x = m'^±/ml , or x = m^o/m 
respectively. The functions F^l'^ are defined as 



A' 



4 



141(1 - 



X 



X) 



122 1 



x) 



(l-x)(l51x2- 335x + 592) 

+6(21x^ - 108x^ - 93x + 50) log 
— 54x(x^ — 2,T — 2) log^ X 

~108x{x^ -2x + 12)Li2(l - x 



8{x'^ - 3x + 2) + (llx^ - 40x + 5) logx 



(37) 



105(1 



-2(x^-2x-2) log^x 

-4(a;2 -2x + 9)Li2(l -x) , (3^ 
(1 - x) ( - 97x2 _ 529x + 2) + 6x2 (I3x ^ j^g^ 



+108x(7x + 4)Li2(l 



X 



(39) 



4(1 



X 



(x + 3) (xlogx + X - 1) + (6x + 2)Li2(l - x) , (40) 



so that they are normalized to unity for x = 1. The logarithms in the first lines of 
Eqs. fl35|) . f p6|) reproduce the result of Ref. |36] for the leading QED-logarithms, 



Aa 



log 



A^SUSY 



times the one- loop result. As can be easily seen, these logarithms are 



negative and reduce the one-loop result by (7. . .9)% for SUSY masses between 
100 . . . 1000 GeV. It is interesting to note that the remaining contributions are 
typically also negative and lead to a further reduction. In particular, for x = 1 
the Ff^-term alone leads to an additional 1.2% reduction of the corresponding 
one-loop contribution. 

For our further discussion of the result we specialize to the case of the MSSM. 
There, the one-loop and the photonic two-loop contributions are given by 



SUSY.IL 



k i,m 



(41) 
(42) 



where the SUSY masses and couplings have to be inserted appropriately, as dis- 
cussed in Sec. [2j 

As a general remark, the photonic two-loop corrections depend on the same 
MSSM parameters as the one-loop contributions, and owing to the structure 
of the analytical results we can expect the overall parameter dependence to be 
similar. In particular, the photonic corrections are proportional to tan /3, just like 



15 



-5 



-5 



= M2/2 
tan p = 50 



(1 = Mj = HoREG = 400 GeV 
Mi = M2/2 
tan p = 50 



-6 



-6 - 




-10 



-10 



400 800 1200 1600 

H = Mj = m|_ ~ = mp ~ = HoFiEQ [GeV] 



2000 



400 



800 1200 



1600 



2000 



(a) 



(b) 



Figure 6: (a) Photonic two-loop corrections, relative to the MSSM one-loop con- 
tributions, as a function of a common SUSY mass scale, see Eq. ( H3!) . (b) The 
same, as a function of mL,fi = mR^p, with fixed M2 = = /xdreg = 400 GeV. 
Our result is shown as a continuous line, the leading-log approximation as a grey 
band. 

the one-loop contributions. Hence, in all our plots we choose a fixed large value 
tan (3 = 50 and plot the ratio of the two- and one-loop contributions. 

Figure [6^ shows the numerical impact of the the photonic corrections, relative 
to the one- loop result, for a simple, generic case|l| We choose the fundamental 
SUSY mass parameters for the Higgsino, wino, left- and right-handed smuon 
equal. 



only the bino mass parameter is determined by the GUT relation Mi = M2/2, 
and tan/3 = 50. The exact result for the ratio of the photonic corrections to the 
one-loop result is denoted by the continuous line. For comparison we also show 
the leading-log result, which, without further knowledge, can only be computed 
with MsusY in a reasonable range. We choose Msusy in the range between the 
minimum and maximum of the mass eigenvalues m^±,o, ?ti^ and represent the 
result by the grey band. 

Fig. [6)d shows a similar plot where M2 = fj, = /ioREG = 400 GeV fixed and 
^L,fi = fnR,fi are varied. 

^In the numerical analysis we follow Refs. (SSHMIE^ and parametrize the one-loop result in 
terms of the muon decay constant G^, i.e. within a^'^®^'^^ we replace na/s^, — > v^G^M^, 
in order to absorb further universal two-loop corrections. This does not influence the explicit 
factor a within a^^^'^^^'^\ which is defined in the on-shell scheme. 



11 = M2= mL,fi = rUR^fi = /iDREG 



(43) 
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Figure 7: (a) Photonic two-loop corrections, relative to the MSSM one-loop con- 
tributions, as a function of /i. (b) Photonic two-loop corrections, relative to 
the MSSM one-loop contributions, as a function of = M2. The other SUSY 
parameters are fixed as shown in the figure. 



Not surprisingly, owing to the additional negative non-logarithmic contribu- 
tions the exact result lies outside the band for the leading-log estimate. The 
leading-log result approximates the exact result best if we choose Msusy as the 
maximum of all SUSY masses (because of the negative sign, this corresponds to 
the lower border of the band). 

In order to understand the behaviour of the MSSM result in more detail 
it is important to note that the MSSM contributions are enhanced by tan/3, 
but this enhancement only affects the terms involving B'"''^; hence the A'"''^- 
terms are comparatively unimportant. The tan /9-enhanced terms can be well 
approximated by mass-insertion diagrams [201 132], with propagating gauginos 
and Higgsinos and insertions of the off-diagonal entries of the chargino/neutralino 
mass matrices. This implies that in the sums over mass eigenstates in f HTl) . f H2l) . 
intricate cancellations take place. The one-loop result depends sensitively on 
mass differences and is mainly determined by the derivatives of the F- \20\ . 

The same discussion can be carried out for the photonic two-loop corrections. 
For the scenario with equal SUSY mass parameters the dominant terms are the 
ones involving and Ff , where (1) = -3/4 and Ff (1) = -45/122. The 
smaller derivative of partially compensates the large coefficient in (13 5p . and 
for this reason the exact result lies only slightly below the leading-log estimate 
for large Msusy- 

Fig. [TK analyses the dependence on the Higgsino mass parameter /i, keeping 
all other SUSY mass parameters fixed. For large or small fi/M2 there is a large 
spread in the chargino mass spectrum, and the uncertainty of the leading-log 



17 




Figure 8: The predictions for for the Snowmass Points and Slopes bench- 

mark scenarios [H]. The wide band corresponds to the present la region of Eq. 
([3]). The narrow band represents the foreseen (SM theory-hmited) improved pre- 
cision if the new g — 2 measurement is carried out given the same central 
value (left) or assuming the deviation vanishes (right). The square points denote 
the full result for a^^^^ including all known two-loop corrections, the diamond 
points are computed without the photonic corrections. 



estimate is large. As the figure shows, in these parameter regions the exact result 
lies in the leading-log band. 

The mass- insert ion diagrams also show that the /i-dependence is non-trivial. 
For small fi, the chargino diagrams dominate. For large n, the diagram with bino 
exchange and left-right smuon transition can dominate [20] — it is the unique 
diagram that increases linearly with /i. But this behaviour is approximately the 
same at the one- and two-loop level, and therefore the ratio shown in the figure is 
almost constant. Nevertheless, for large /i positive photonic contributions start 
to partially cancel the leading logarithms. 

Fig. [7b shows the behaviour if all chargino and neutralino masses are varied 
together, fi = M2 and Mi = M2/2, with fixed smuon/sneutrino mass parameters. 
Similarly to Fig. [7^, positive two-loop contributions can become important at 
large /i, and the exact result can even lie above the leading-log band. 

Finally we consider the results for the SPS benchmark points [H]. E.g. SPSla 
and SPSlb lead to a^^^"^ close to the observed deviation Q. SPSla has tan/3 = 
10 and quite small SUSY masses, SPSlb has tan/3 = 30 and slightly larger SUSY 
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masses. The results for these points are, in units of 10 





„SUSY,1L 


„SUSY,2L leading log 


^SUSY,2L (7) 


SPSla 


30.49 


-1.93. . . - 2.32 


-2.18 


SPSlb 


33.34 


-2.27. . . - 2.63 


-2.59 



Fig. E] shows a graphical distribution of the ten SPS benchmark predictions 
for a^^^ , computed with all known one- and two-loop corrections, and compared 
with the experimental and SM value. The discriminating power of the current 
and an improved determination is evident. The plot also displays the results 
obtained neglecting the photonic corrections. For several benchmark points, the 
difference is as large as one sigma of the future Aa^(exp — SM). 



6 Matching the result to a full MSSM calcula- 
tion 

Our results have been presented in such a way that they can be easily utilized 
as building blocks in a full calculation of all two-loop corrections to aSUSY,iL 
the MSSM. This is useful because two technical problems are concentrated in the 
photonic corrections: the appearance of IR divergences and of terms which are 
not suppressed by powers of divided by a heavy mass of the order of the weak 
scale or SUSY scale. 

As discussed in Sees. [H 121 the SUSY two-loop corrections to SM one-loop 
diagrams are already known. What remains to be computed are the two-loop 
corrections to SUSY one-loop diagrams, denoted by a^^^'"^^^^^ in Ref. [20]. These 
contributions are defined as the two-loop diagrams where the /z-lepton number is 
carried by a /i and/or z/^ line, plus corresponding counterterm diagrams. 

Similar to our photonic corrections, these full contributions can be decom- 
posed as '^''^s^"" + a;^*'Q™^^" + a;^*'SUSYfuii ^ ^ct,rem_ gg^g ^2Ldiagfuii ^^^^^^ ^he 

genuine two-loop diagrams and the counterterm contributions have been split 
into the QED counterterm diagrams of the form in Fig. HI the SUSY counterterm 
diagrams of the form in Fig. O and all other remaining counterterm diagrams. 
Hence we can write 

^SUSY,2L(b) _ SUSY,2L(7) 

_^^^2Ldiagfull_^^2Ldiag(7)^ 
_|_ ^^ct.QEDfuU _ ^ct.QEP-^ 
+ (aJ'SUSYfull_^^ct,SUSY) 

+ (aj'--), (44) 
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where the appropriate summation over the chargino, neutrahno and smuon in- 
dices of each term is imphed. 

Setting up the full calculation in this way has several technical advantages. 

• In the difference of the genuine loop diagrams (2nd line of Eq. dS])), simply 
all the photonic diagrams of Figs. [21 |3]drop out. In other words, only the 
non-photonic two-loop diagrams need to be evaluated. As an advantage, 
these are all individually infrared finite and suppressed by divided by a 
heavy mass squared. 

• The difference of the QED counterterm diagrams (3rd line of Eq. (I44l) ) van- 
ishes. The reason is that the appearing QED counterterm insertions and 
renormalization constants, see Eq. fl22|) . have to be defined in the same way, 
in the on-shell scheme, in the full theory and for the photonic corrections. 
Therefore, the QED counterterm diagrams, which are again partially in- 
frared divergent, need not be re-evaluated in the full calculation. All coun- 
terterm diagrams that remain to be calculated are individually infrared 
finite and suppressed by divided by a heavy mass squared. 

The only subtlety arises in the 4th line of Eq. (jBl), in connection with the 
SUSY counterterm diagrams of Fig. [5l The renormalization constants ( l23l) ap- 
pearing within these diagrams are different in the full MSSM and in our cal- 
culation. In the full theory, must also be defined in the on-shell scheme, 
but further diagrams contribute to it. Mass and field renormalization constants 
for the neutral particles do not vanish any more. Finally, for the SUSY masses 
and couplings not even the same renormalization scheme can be used in the full 
theorylfl 

Nevertheless, it is possible to choose any desired MSSM renormalization 
scheme and to compute the full MSSM counterterm diagrams a^t.susYfuii tY^oX, 
scheme. Then simply our counterterm result a^*'^^^^, which has been given in 
Eq. (127]) . has to be explicitly subtracted. The choice of the regularization scheme 
DREG or DRED for the counterterms has to match the choice for the two-loop 
diagrams. In this way, the full result corresponding to the desired MSSM renor- 
malization scheme is obtained. As discussed in Sec. HI a^t'SUSYfuU ^jjj contain the 
same large QED logarithm arising from (5m^ within the SUSY coupling renor- 
malization that is already contained in a^*'^^^"^. Therefore, the difference in the 
4th line of Eq. f l44p is free of large QED logarithms. 



^In the MSSM, supersymmetry together with SU(2)xU(l) gauge invariance imphes correla- 
tions between the SUSY masses and the different couplings c, n, which must be reflected in the 



renormalization scheme, see e.g. |49H51) . But they cannot be reflected in the purely photonic 
corrections since these are not invariant under the full symmetry of the MSSM. 
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7 Conclusions 



In the present paper the photonic two-loop corrections to the muon magnetic 
moment in the MSSM and a wider class of models have been evaluated exactly. 
The photonic corrections are defined as the two-loop diagrams which contain 
a photon loop attached to a SUSY one-loop diagram, plus the corresponding 
counterterm diagrams. The counterterms are defined in the "on-shell muon mass 
scheme" , which is natural for the models considered here and in Ref. [3S] — all 
QED quantities and all masses are renormalized in the on-shell scheme, but the 
purely high-scale parameters are MS-renormalized. 

Our result reproduces the large logarithms of Ref. [36] and provides the exact 
result for the additional subleading logarithms, dilogarithms and non-logarithmic 
terms. The leading logarithm has an intrinsic uncertainty because it could be 
evaluated with a small, large, or intermediate SUSY mass. For the typical SUSY 
scenarios considered here, the leading logarithm amounts to around —7% of 
a^^^^'^^ if the smallest SUSY mass is used; the additional terms are in the range 
(0.5 ... — 2)% and thus typically lead to a further reduction of a^^^^'^^. How- 
ever, for large /x, the additional terms can have a positive sign and can partially 
compensate the leading logarithm. 

It is interesting to compare the photonic two-loop corrections to other known 
SUSY two-loop contributions to the muon magnetic moment. Another universal 
two-loop correction is the (tan/5)^-correction arising from a shift of the muon 
Yukawa coupling [37]. In a large part of the MSSM parameter space, particu- 
larly for approximately degenerate SUSY masses, the photonic and the (tan/3)^- 
corrections are the largest two-loop effects. While the photonic corrections are 
negative, the (tan /3) ^-corrections are positive (for positive a^^^"^'^^) and can over- 
compensate the photonic corrections for large tan (3. The SUSY two-loop correc- 
tions to Standard Model one-loop diagrams [331 [M] amount to around 2% of 
a^^^^'^^ for degenerate masses, but in special cases with large mass splittings 
or in cases where the one-loop contributions are suppressed, they can become 
dominant. 

Importantly, with the exact two-loop computation the theory error arising 
from unknown photonic corrections has been reduced. The remaining theory error 
due to unknown photonic three-loop corrections can be estimated by comparing 
with the electroweak contributions in the Standard Model. There, the photonic 
three-loop corrections amount to only 1% of the photonic two-loop corrections. 
Based on that, we can estimate the unknown photonic three- loop corrections to 
qSUSY,il Yqss than 0(0.1 x 10"^°) and thus negligible. 

Nevertheless, in order to make full use of the expected and intended improve- 
ments of a^^P and for SUSY phenomenology, the SUSY theory error should 
be further reduced. Our calculation can also help in computing the remaining 
SUSY two- loop contributions to a^. If the full calculation is organized as in Eq. 
(jH]), only non-photonic diagrams and non-QED counterterm diagrams remain to 
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be computed. These are all infrared finite and suppressed by the required 
factor and contain no further large QED-logarithm. 
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